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Nomenclature 


Cx. Cy, C z 

ce 


e, e x 

f, fy 

F x 


H 


K, 


1J 


Fluid velocity components along x y z 
same as Cy 

Fluid Damping Matrix for turbine displacements 

Turbine eccentricity 
Tangential force per unit length 
Net force on turbine in direction of offset 
Net force on turbine perpendicular to offset 

Blade height 

Huid stiffness matrix for turbine displacements 
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Structural stiffness for turbine displacements 
Rotor mass 

Pressure, total pressure 
Turbine torque 
Turbine mean radius 
Turbine wheel speed (wR) 

Axial coordinate 
Tangential coordinate 
Radial coordinate 

Force coefficient along offset (Eq. 8) 

Force coefficient perpendicular to offset (Eq. 8) 

Sensitivity of tangential force to relative gap, (J — — 
Sensitivity of local efficiency to relative gap, p = - 


3f„ 


3(5 /H) 
3fi 


3(5 /H) 

Blade tip gap 

Row coefficient, C x /coR 

Phase angles for fluid force (Eq. 10) and wall pressure (Eq. 1 3) 

Ruid density 
Azimuth angle 

Angular frequency of turbine spin 
Whirl angular frequency 

Rotor damping factor 


1. Introduction 

Th . pvietinff laree body of literature on turbomachine tip leakage has been 

promote forward whirl, leading to rotordynamic instability. 

Both Alford and Thomas showed that, if the tangential force f per unit length is 
assumed to vary linearly with the ratio 5 / H of local tip gap to blade height 




( 1 ) 


and if the shaft is offset instantaneously by e x along the o x direction, then a cross- 
force F y equivalent to 


F 

y 2 R H 
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arises where Q is the turbine torque and R its mean radius. The factor P has 
become known in the U.S. as the Alford coefficient which is quoted as varying between 0 
and 6. The German literature uses a factor called the "excitation coefficient", k 2 , 

which is equivalent to P/2 . Implicit in the Alfonl/Thomas model is the assumption 
that the flow remains perfectly uniform up to the offset turbine, as well as downstream of 
it, so that only the local efficiency, as determined by the local tip gap, is of concern. Two 
of the consequences of these models are (a) the absence of a direct force component (in 
the offset direction, and (b) the absence of frequency dependence (no damping). 

If only a cross-coupled coefficient (K X y) is generated, it can be easily shown that 
the damping factor required for stabilization against it would be 

£ . =ii^i=P_9_ m 

mul 2 K c 4 RH K„ (3) 

This can be very substantial in high-pressure turbines. For the SSME H-? 
turbopump, Q = 12800 Nm, R = 0.129m, H = 0.023m, and Kq = Mrot 

Q„ = 1.9xl0 2 * * * * * 8 Nm, giving £ min = 0.019, or a log decrement of 12%. 

Relatively little work on these forces has been done since the pioneering efforts 
of Alford and Thomas. Thomas' collaborators at the T.U. Munich produced the most 
detailed experimental data. Urlichs (3) used a relatively low power facility, with blade 
Reynolds numbers below 10^, and measured cross-forces mainly on shrouded turbines, 
although one unshrouded case was also tested. He identified the shroud seal as the major 
contributor to the cross- forces. The mechanics of these seal-related forces has been more 
clearly elucidated since (see for instance Refs. (4), (5), (6)), and is distinct from the 
uneven work extraction of Alford-Thomas effect. In his unshrouded tests, Urlichs noted 
a cross-force reduction with increasing mean tip gap, and an increase with axial stator- 
rotor spacing. The measured forces were roughly compatible with the simple Alford 
argument. Wohlrab (7), used a larger, pressurized air turbine, capable of stator leaving 
Reynolds numbers up to 5x10^, but tested only shrouded turbines. Strong non-linearity 
of force vs. displacement was noted in the cases with smaller forces, which raises 
questions about accuracy. Limited dynamic testing was accomplished as well. As in 
Ref. (3), the main mechanism in these tests was through seal pressure effects, rather than 
through uneven work extraction. Vance and Laudadio (8), did some very low power tests 
on a fan with statically moveable casing, and reported measured cross-forces, but no 
aerodynamic data. Ehrich (9) has recently inferred Alford forces from compressor 
efficiency test data, and argues that these forces become backward-whirling at pressure 
ratios above the normal operating point. 

We have performed extensive force measurements on an unshrouded turbine 
identical to the first stage of the Shuttle LH turbopump, and also on a shrouded derivative 
of it. These were supplemented by flow field measurements and theoretical analysis to 
clarify mechanisms. In this paper, we will describe the test facility (Sec. 2), and present 
the force and flow data (Sec. 3-6). 


2. Experiment al Apparatus 

The test facility is a pressurized closed loop, filled with Freon 12 gas, and 

equipped with a gas blower, heat exchanger/cooler, removable test section, power 

extraction generator and data acquisition system. Nominal operating conditions are 2.0 

atm mean pressure and 4.5 Kg/sec flow rate. Flow rate is mainly controlled by a manual 

series valve, with help from a bypass valve at low flow. Speed is controlled through 

generator excitation control. 
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The test section area is shown in Fig. 1. The upper section (12) of the casing can 
be rotated and carries the stator and the hub, as well as a variety of flow probes. The 
turbine shaft connects to the turbine via a four-post rotating dynamometer (14), and is 
supported by two bearings. The bearings are carried in a heavy cylindrical structure 
which can be translated sideways by means of four stiff rods, two on each side. Static 
offsets are achieved (to an accuracy of ± 0.5 mil) by insertion or removal of calibrated 
shims (11). The dynamometer can sense all components of force and torque on the 
turbine, and the signals from its 9 strain gauge bridges are carried through the shaft to the 
slip-ring assembly (23). These signals were sampled 32 or 64 times per revolution, on a 
pattern which was phase-locked to the rotor by means of signals from an auxiliary 
encoder (18). The data were ensemble-averaged over 128 or 256 revolutions, to reduce 
low-frequency noise, and were then numerically projected on fixed axes to extract the DC 
forces of interest (F x along the displacement axis, Fy perpendicular to it). Despite 
sizable second harmonic contamination from the flex joints in the shaft, these DC 
components were extracted with an accuracy and repeatability of ±0.05 /b f (the forces 
F x and Fy ranged to a fewfb f ). 

The characteristics and design parameters of the unshrouded test turbine arc 
summarized in Table 1. Due to the low pressure ratio, compressibility effects are 
minimal, and similarity would be ensured by matching flow and work coefficients, and 
Reynolds number. The first two are indeed matched to the SSME LH turbine. The 
Reynolds number (based on stator exit velocity and blade height) is 1.4xl0 6 , compared to 
5.6x10^ in the SSME turbine; since both are well above transition (10^ -2x10^) no 
significant difference is expected. A limited test of this insensitivity was provided by 
comparison runs at 1 atm and 2 atm loop pressure, in which no difference was detected in 
the nondimensional cross-force characteristics. 

After the unshrouded turbine tests, the same turbine was modified by removing 
the outer 6.49 mm (out of a blade height of 22.9 mm) and installing a continuous shroud 
band with a 2-ridge labyrinth seal (as shown in Fig. 2). Due to the partial flow blockage 
of the unrecessed seal, the optimum flow rate at the design speed of 3440 lpm dropped to 
3.16 Kg/sec, the optimum efficiency dropped to 74%, and the corresponding pressure 
ratio was reduced to 1.14. 

The test matrix is summarized in Table 2, with the geometrical notation 
contained in Fig. 3. 


3. Force Data for the Unshroud ed Turbine 

Typical force vs. displacement plots are shown in Fig. 4, which corresponds to 
Configuration 1 of Table 2, at its design condition. The plots for Configurations 2-5 are 
qualitatively similar. In these graphs, the negative F x slope indicates a restoring direct 
force, while the positive Fy slope indicates a forward-whirling cross-force. Each test was 
repeated three times, and all test results are shown to illustrate the degree of repeatability 
of the data. The fact that the forces are not exactly zero at zero eccentricity is due to a 
combination of casing out-of-roundness and positioning error. Despite the relatively 
large offsets (± 15 mil on a mean gap of 27 mil), the Fy data show no departure from 
linearity. By contrast, the F x data show in all cases a slight s-shaped curvature, with the 
slope increasing with eccentricity. 

As noted in the Introduction, the simple Alford theory would predict F x = 0. 
However, the data show IF X I of the same order as IFyl, and a different mechanism, ora 
variation on Alford's postulated mechanism, must be involved. 
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The results of the measurements for all configurations are reported in 
nondimensional form in Table 3. The coefficients ax and ay are obvious 
generalizations of Alford's b: 


__2F i R_ . 2F y R 

~~ Q(e / H) ’ “ y Q(e / H) 


(4) 


where the notation has been changed from b'j which is strictly the sensitivity of blade 
tangential forces to tip gap, to a, which is a measure of the cross-forces, and may or may 
not be equal to b (indeed, there is no b counterpart to a x ). The forces at e=0 have been 
subtracted from F x and Fy in calculating a x and ay. 

Table 3 indicates a general increase of a y with speed at a fixed flow rate, or a 
decrease with the flow coefficient <|> = C x /(wR). This is displayed in Fig. 5. The 
theoretical curve shown is from an xz actuator disk theory which will not be discussed in 
this paper (see Ref. 10). 


A second trend in the data is a substantial increase of the force coefficient (both 
la x l and ay) as the mean radial gap is reduced. This can be seen by comparing 
Configurations 2 (gap 3% of blade height) and 4 (gap 1.9%), both with the widest axial 
hub gap d', and also by comparison of Configurations 3 and 5, similarly related, but with 
a narrower axial hub gap. Averaging over the various speeds, the effect amounts to a 0.6 

decrease in a y (19%) per 1% increase in 5 / H (wide axial gap), or a 0.7 (21%) per 1% 
for the narrow axial gap. Of course, since only two gap values were tested, there is no 
confirmation of the linearity of this effect. The theory (Ref. 10) predicts a much smaller 

effect of 6 on ay, and there are some indications from the flow data (Sec. 4) that viscous 
flow effects in the narrow tip flow passages may be responsible. This trend had been 
previously reported by Urlichs (3) as well. 


Configurations 2 and 3 differ only in the hub axial gap d', while Configuration 1 
differs from 2 and 3 in both, d', and the stator-to-rotor blade spacing, d. All three 

configurations have 5 / H = 3%. Thus, if we postulate a linear variation of ay with d 
and d’, the data in Table 3 can be used to extract the separate sensitivities of a y to these 
gap values. Similarly, comparison of Configurations 4 and 5 can yield the d' sensitivity 
for the cases with a narrow radial tip gap. The results (Table 4) are inconclusive for d, 
but are unambiguous as to sign and general magnitude for the hub gap d'. We have not 
been able so far to find a satisfactory explanation for this effect. Opening the gap d' 
should have the direct effect of reducing the pressure non-uniformity in the stator-rotor 
space, and to the extent that this nonuniformity contributes to the cross-force (see Sec. 4), 
this would indeed reduce ay. On the other hand, these pressure nonuniformities also 
redistribute the upstream flow in a manner which tends to dampen the Alford effect, as 
our x-y actuator disk analysis (Ref. 10) makes clear. The net result must then depend on 
the balance of these two effects. A more complete analysis of the flow data, and 
additional theoretical development are needed in this area. Uriichs (3) found the opposite 
trend (ay increasing with axial gap), but his geometry was such that belli d and d' were 
varied simultaneously. 


In addition to the direct force measurements, both upstream and downstream 
flow fields were explored through a set of static wall pressure taps and directional probes. 
More details about this instrumentation and flow data can be found in Ref. 10. The 
surveys upstream of the stator indicated very small departures from tangential uniformity. 
On the other hand, the wall pressure surveys downstream of the stator do show a well- 
resolved nonuniformity. Fig. 6 shows the pressure pattern for the inter-blade row region 
and the rotor blade tip region. For the region between stator and rotor a pressure 

fluctuation amplitude of about 0.0027P to = 0.22 pC* 0 / 2 = 0.028p(coR) 2 for e/H = 
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0.019. The pressure minimum is about 25* ahead of the maximum gap location The 
pressure nonumformity increases in amplitude as one moves downstream over the rotor 
blades, particidarly in their narrow-gap sectors, but it then returns downstream to the 
same level as between rotor and stator. Thus, the turbine eccentricity is felt upstream and 
downstream of it, through potential effects, with a length scale of the radius R, and, from 
limited probing data, spanning the passage depth as well. 

At 3 chord lengths downstream of the turbine, the velocity data were obtained to 
a depth of 75% span, using a 3 hole probe. The tangential distribution is shown in Figure 
7, where the mean (centered turbine) tangential velocity has been subtracted out. Based 
on consistency among various probes and also with the dynamometer data, the errors 
appear to be under io%. The positive tangential velocity values seen near the bigger gap 
region (180 ) indicate undertuming by die rotor blades. The magnitude seems to be 
uniform in the outer 5% of the span, gradually decreasing towards the hub. 

Thus, two main sources of lateral force on the turbine due to the flow properties 
can be identified. 1) a tangentially non-uniform flow turning which leads to uneven work 
extraction; and 2) a non-uniform static pressure distribution acting on the turbine hub. 
Both forces can be integrated around the perimeter to a direct force and a cross force! 
Integration of the pressure nonuniformity is straightforward. We assumed the pressure 
acting on the hub is that observed before and after ther rotor, since the excess 
nonuniformity over the rotor blades appears to be a near-tip effect. For the work-defect 
forces, we used the Euler turbine equation locally to calculate an azimuth-dependent 
tangential force per unit length 

f y =Jo I pC x AC y dr (5) 

where C x is taken to be the mean axial velocity, and the flow turning ACy is from our 
velocity surveys. 

The results for Configurations 4 and 5 are shown in Table 5 , where we report the 
separate pressure and work-defect contributions to ax and ay, their sum, and for 
reference, the dynamometer-measured data. Similar results were also obtained (for all 
configurations) using velocity data at 1.5 chords downstream of the rotor (Ref. 10). Most 
of the direct force is due to the pressure force while the cross force is caused by both the 
pressure force and the blade force. The cross force coefficients, a y , are larger than 
expected for this low-reaction turbine on the basis of work-defect arguments. Clearly, 
the pressure component, not previously identified, is an important factor. 

All of the data reported here are for static offsets. An attempt at measuring cross 
forces under dynamic conditions (using an inertial shaker) was unsuccessful due to 
excessive vibratory noise. The simple theory of Alford predicts zero damping, while a 
more complete model (Ref. 10) indicates significant damping due to lag in tangential 
flow redistribution. Given the large magnitude of the statically measured forces, 
measurements under dynamic conditions should be undertaken in order to complete our 
understanding of Alford forces. 


4. 


The turbine was modified, as described in Section 2 of this paper, for the 
shrouded case, and dynamometer and flow data were obtained. As in the unshrouded 
cases, the forces scale linearly with eccentricity, and a sinusoidal pressure pattern with a 
large amplitude develops over the shroud band. Table 6 shows the excitation coefficients 
from the dynamometer data, and the pressure data at the design speed. Compared to the 
unshrouded cases, the excitation coefficients are larger by a factor ranging from 1.5 to 
2 . 0 . 
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The nonuniform pressure distribution produces both direct and cross forces 
which are smaller than those measured with the dynamometer. For this configuration, no 
measurements were taken of the velocity perturbations, because, due to the presence of 
the unrecessed shroud, the flow was very turbulent downstream of the rotor. It it is 
assumed that the work defect mechanism contributes pnmanly to the cross force as in the 
unshrouded cases, the discrepancy in the direct force excitation coefficient cannot be 
explained. Furthermore, the small difference in the cross force excitation cc^fficients for 
the shrouded case suggests that the pressure effect, instead of the work defect effect, is 
primarily responsible for the cross force. Lastly, the pressure non-uniformity, which can 
be detected upstream and downstream of the shroud, shows again a flow redistribution on 
the scale of the turbine radius. 

A linear labyrinth seal model based on the work of Millsaps (6) was extended to 
include the effects of non-uniformities in the flow both upstre^ and downstream of the 
shroud and was used to analyze the shrouded data. It was found that the non-^formihes 
have a large effect on the model's predictions, essentially doubling the magnitude : of both 
the direct and the cross force, but even after this correction, the model 4 nd erpredicts die 
forces by about 40%. No complete theory exists of a seal interacting with the flow field 
of the turbine blading. Palczynshki (11) discusses, in greater detail, the work on the 
shrouded turbine. 


g. Conclusions 


This work has confirmed the existence of the destabilizing forces suggested by 
Thomas and Alford. The general scaling and order of magnitude are also consistent with 
their insights. However, some new effects in the unshrouded cases include the following. 


1) In addition to a non-uniform work extraction, a non-uniform pressure distribution 
also exists. 


2) The forces increase significantly as the mean tip gap is reduced, confirming an earlier 
result. 


The lateral forces in shrouded turbines can be larger than those in unshrouded 
turbines. This is due to a large nonuniformity in the seal gland pressure, which now 
dominates over the work defect contribution. _ 

Alford force measurements under dynamic conditions are recommended in order to 
investigate the possible existence and magnitude of a damping component of this torce. 
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How coefficient, <p 


Work coefficient, 

Stator exit angle 
Relative rotor inlet angle 
Rotor exit angle 
Absolute exit angle 
Degree of reaction 
Rotor mean radius, cm (in) 
Number of rotor blades 
Rotor blade height, cm (in) 
Rotor blade chord, cm (in) 
Rotation Rate, rpm 
Axial flow velocity, m/s (in/s) 
Mass flow rate, kg/s (slug/s) 
Inlet pressure, kPa (psi) 
Inlet temperature 
Pressure ratio 
Efficiency 



_ Table 1 

Design Parameters for SSME Fuel Turbopump Fi ret Stage 
and the Test Facility Turbine 


uel Turbopump Rrst 
Stage 




70* 

43.9° 

60° 

-3.1* 

0.216 

12.88 (5.07) 
63 

2.17(0.854) 
2.21 (0.870) 
34,560 
262 (10,300) 
71.8(4.92) 
34,950 (5069) 
1053 (1436) 
1.192 
0.821 




-60° 

-3° 

0.216 

12.88(5.07) 

63 

2.17(0.854) 
2.21 (0.870) 
3440 
26 (1020) 
4.48 (0.307) 
224 (32.43) 
300 (80) 
1.231 
0.80 
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Table 2 The Test Matrix 


w/wjcs 

(%) 
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* Table 4 

Sensitivity of Cross-Force to Axial Stator-Rotor Spacing (d) 

(See Figure 6 for definitions) 


and Hub Axial Gap (d') 


S/H 


da y / d(d/c ) 

Ba y / d(cf /c) 

0.030 

0.030 

0.030 

r o.7 i 
1.0 
1.1 

— m — 

-1.21 

-1.00 

T75 

-0.44 

-1.44 

0.019 

0.019 

0.019 

077 “i 
1.0 
1.1 

• i i 

332 — 

-3.13 
- 2.37 


Table 5 Contnbutions from work defect (wd) and pressure (p) to the force coefficients 
a x, cty at design conditions, and comparison to dynamometer data. 


Configuration 

3 


( a x)wd 

”TT" 




5 


-.6 


-2.3 


( a x)total 

-3.2 

-2.9 


(«x)dyn 

~TW~ 

-3.47 


ter and 
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(1) Flow straightener 

(2) Screen 

(3) Main loop piping 
^4) Range 

(5) Rcxible insert 

(6) Liner 

(7) Snubber bearing 

(8) Snubber support 

(9) Teat turbine 

(10) Row-smoothing shield 

(11) Shims 

(12) Rotatable casing 

(13) Stator blades 

(14) Rotating dynamometer 

(15) Bolts to secure shaft 

(16) Turbine -offsetting rods 

(17) Upper flex joint 

(18) Optical encoder 

(19) Intermediate shaft 

(20) Double-acting seal 

(21) Flexible insen 

(22) Pivoting bearing 

(23) Slip ring assembly 

(24) Lower flex joint 


Figure 1: Schematic of the turbine test section. 




Figure 3: Schematic of turbine's major dimensions of 
interest. 


Figure 2: Shroud-seal for turbine. Dimensions in 
mm. 
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